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ABSTRACT 

The  research  performed  under  Grant  AF-EOAR-63-76, 
is  reported  in  three  main  parts. 

Part  1 

This  mainly  covers  the  material  in  the  original 
proposal.  A  rough  model  of  the  rocket  combustion  process 
demonstrates  that  if  a  modulated  flow  of  fuel  could  be  produced  at 
the  injector  plate,  a  repetitive  and  corresponding  variation  of 
combustion  pressure  would  result,  and  form  a  driving  mechanism  for 
high  frequency  combustion  instability.  It  must  be  noted  that  at  the 
frequencies  involved,  the  inertia  of  the  fuel  feed  lines  would  hinder 
the  total  fluid  column  behind  the  injectors  from  following.  Thus 
the  fuel  flow  modulation  must  come  from  a  modulation  of  the  discharge 
velocity  from  each  orifice,  for  example,  by  a  change  in  the  pressure 
at  the  plane  of  the  orifice  exit  (this  pressure  determines  the  dis¬ 
charge  velocity  for  a  given  fuel  and  orifice). 

The  high  *  g*  vibrations  to  which  rocket  engines  are 
subjected  both  on  the  test  stand  and  in  flight,  suggested  that 
transverse  mechanical  vibration  of  the  injector  plate  could  induce 
this  fuel  flow  modulation. 

Part  2 

Initial  work  on  the  effect  of  low  *  g’  vibration  on 
injector  behaviour  was  reported  at  the  l6th  AFOSR  Meeting  on 
Nantucket  Island  and  in  the  June  1964  edition  of  the  British 
Institute  of  Applied  Physics. 

It  has  been  established  that  vibration  induces  velocity 
modulation  in  the  jet  emerging  frcm  the  injector.  At  low  *  g’  the 
minute  modulation  merely  serves  to  trigger  the  capillary  type 
instability.  At  high  *  g’  (greater  than  100),  the  velocity  mod¬ 
ulation  is  sufficient  to  result  in  a  predominant  bunching  effect  in 
the  liquid  jet.  Spectacular  disc  formation  on  the  jet  is  the  result. 
Discs  on  jets  frcm  neighbouring  holes  can  be  in  phase  or  out  of  phase 
depending  on  their  alignment  with  respect  to  the  line  of  vibration. 

’Bunching*  occurs  -whether  the  jet  is  laminar  or 
turbulent,  low  velocity  or  high  velocity. 

Rayleigh*  s  original  analysis  has  been  modified 

(a)  to  specifically  incorporate  velocity  modulation  as  a  trigger 
mechanism,  and  to  cover  second  order  effects. 

(b)  to  include  the  ’bunching*  effect  at  least  in  the  transition 
region  where  it  has  not  become  completely  dominant  over  the 
capillary  effect. 
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then  to  couple  with  an  acoustic  mode  of  the  chamber  to  become 
self  sustaining. 

A  tangential  mode  of  vibration,  by  including 
space-wise  phase  shift  in  the  jet  modulation  across  the  injector 
plate,  can  initiate  the  tangential  mode  of  instability.  Thus, 
modes  of  •ombustion  instability  can  be  related  directly  to  modes 
of  vibration. 

Part  3 


The  effect  of  the  tangential  mode  of  instability 
on  the  boundary  layer  heat  transfer  coefficient  is  considered. 

The  wall  vibrations  set  up,  cause  driven  vortices 
of  the  Taylor-Goertler  type  with  axes  parallel  to  the  cylindrical 
chamber  axis. 

The  increase  in  the  heat  transfer  coefficient  with 
respect  to  that  in  the  normal  laminar  layer,  i3  derived.  A  small 
effect  was  found. 

The  presence,  however,  of  rotational  motion  of  the 
chamber  gas,  associated  with  the  tangential  mode,  leads  to  an 
unsteady  vortex  motion  in  the  boundary  layer,  which  builds  up 
rapidly.  These  vortices  resemble  those  observed  to  be  formed 
between  two  concentric  rotating  cylinders. 

This  vortex  growth  causes  a  very  significant  increase 
in  the  heat  transfer  coefficient  which  occurs  in  a  very  small  time 
(less  than  10  sec) .  Upper  and  lower  limits  of  the  effect  are 
established. 

The  increase  in  heat  transfer  coefficient,  combined 
with  the  increase  in  rate  of  heat  release  associated  with  the 
tangential  mode  in  a  plane  close  to  the  injectors,  will  normally 
cause  rapid  chamber  wall  burn-out. 

It  is  considered  that  turbulence  in  the  rocket  engine 
boundary  layer  is  unlikely,  and  that  separation  is  not  present. 
That  is  to  say,  a  well-ordered  and  effectively  laminar  layer 
exists,  even  though  vortex  formation  has  occurred. 

Vorticity  induced  in  the  boundary  layer  due  to  shear 
flow  in  the  fluid  external  to  it,  is  known  to  have  a  small  effect 
on  the  heat  transfer  coefficient  (certainly  much  less  than  that 
due  to  full  vortex  formation  the  boundary  layer). 

Menu 3,  the  only  source  of  appreciable  increase  in 
the  coefficient  is  that  due  to  the  rapid  build-up  of  the  Taylor- 
Goertler  vortex  formation  with  free  stream  rotational  motion. 

The  work  points  to  a  possible  way  of  preventing 
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”A  driving  Kechanism  for  High 
Frequency  Combustion  Instability 
in  Liquid  Fuel  nockot  Engines" . 


Section  l .  Introduction 


'/■»«.  arable**'  of  combustion  pr-.ssuro  oscillation  in  li  -uitl- 
fu-I  rocket  a.  tor  ope  rrLi.cn  L.-e  lor.''  n.e.-i  L-.o  subject  of  theoretical 
and  experimental  investigations. 

She  low  frequency  (leas  Urn  2CC  c.n.o.)  tyoo  of  oscillation, 
known  ns  " chuggii i-f j* ,  has  ixjcn  t,  .orcug'nly  analysed  and  the  problem 
solved  -  sou  Jrocco,  5th  Combustion  dym xosium,  p.164. 

This  yipor  is  concerned  with  the  more  complex  (and  more 
destructive)  high  frequency  oscillations,  covering  a  range  X'rov 
about  1,000  to  6,0C0  c.p,3.  .juch  oscillations  can  resonate  with 
the  acoust'  cal  varies  of  too  combustion  chamber.  Longitudinal, 
tangential  and  radial  oscillating  modes  have  been  ooservod. 

Pick  ford  and  Peoples  (1)  have  put  forward  a  very  cum- 
prehensive  treatment  of  the  inherent  stability  of  the  combustion 
orocess.  however  their  work  involves  the  postulation  of  <a  basic 
pressure  perturbation  which  i3  not  define.!.  This  results  in  the 
use  of  a  rather  incomplete  oerturbntion  index. 

That  injection  is  an  innortant  factor  in  high  intensity 
combustion  systems,  whereby  the  -itemization  process  is  used  to 
create  the  maximum  practical  surface  area  of  liquid  fuel  to  aid 
rapid  evaporation,  na3  been  recognised  for  many  years. 

Under  conditions  of  high  temperature  and  pressure, 
chemical  kinetics  heretofore  was  considered  as  not  being  a  con¬ 
tributing  factor.  Characteristic  reaction  times  are  of  the  order  of 
1CT7  sec,  in  ccmp.-irison  with  character:  otic  chamber  frequencies  with 

o 

periods  in  the  region  of  10  sec.  Lost  of  the  residence  time  in 
the  combustion  chamber  is  saent  in  the  physical  processes  of  liquid 
atomization,  drop  evaporation  and  gas  ohase  nixing,,  prior  to 
reaction,  all  of  these  factors  are  closely  related  to  the  injector 
design , 

Several  theories  of  combustion  system  performance  have  been 
oostulated  on  the  basis  of  ovaooration  as  the  only  ruto  con  troll  in:* 
process.  However,  practical  studies  indicate  differences  between 

the  predicted  system  behavious  and  that  actually  observed . 

In  any  case,  it  does  not  appear  tiv*t  the  evaporation  process  alone 
can  explain  the  role  of  the  injector  as  an  initiating  instability 
mechanism. 

Progress  in  theory  and  experiment  then,  so  far  has  only  led 
to  certain  rough  rules  /nich  act  as  guides  in  the  desigp  stage  of 
rocket  engines.  Pressure  oscillations  can,  and  still  do,  occur 
at  random  io  even  the  moot  conservatively  .L signed  units,  lending 
at  least  to  severe  impairment  of  o aeration. 

The  establishment  of  the  basic  cause  of  f  .13  phenomenon 
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Tho  nrescncc  of  random  combustion  nressure  fluctuations  in 
a  smoothly  running  engine  is  easily  demonstrated •  Tue  aml.itud  xs  are 
s olden  greater  than  about  one  or  two  percent  of  Lhc  steady  pressure. 
These  fluctuations  could  be  a  consequence  of  discontinuous  energy 
injection  Into  the  system.  Any  mechanism  which  would  synchronize 
events  at  all  injectors,  could  t  en  cause  considerable  pressure 
oscillations . 

CD 

It  is  known  that  both  on  the  te3t  stand  and  in  actual 

flight/5^  tho  rocket  engine  can  be  subjected  to  hi$i  1  g*  higjh 

frequency  vibration  (dee  loacadix  1). 

In  tho  140.CC0  lbs  thrust  Atlas  type  booster,  taken  later 

(6) 

as  an  example,  theory  predicts  a  30  1  g* ,  10CC  c.p.s.  response  to 
the  conponont  of  tho  engine  noise  field  arriving  at  the  vcuiclo 
surface.  On  test  stands  much  higher  * g*  values  have  been  observed. 

It  is  contended  that  this  vibration  will  cause  periodic 
(at  the  applied  frequency)  variation  of  the  fuel  flow  rate  from  an 
injector.  This  “hunching"  of  fluid  moreover,  ‘./ill  be  3ync  .ronized 
at  each  injector  (all  the  holes  will  oe  synenronizod  provided  tho 
injector  plate  is  vibrating  in  tho  n  -  1  mode) .  •Pcperimcntal 
verification  and.  investigation  of  t-is  phenomena  has  occupied  most 
of  the  work  during  this  grant  and  will  bo  rcoortcd  in  .’art  2. 

Thus,  with  synchronization,  the  fuel  flow  rate  entering  the 
combustion  zone  will  be  modulated  at  the  vibration  frequency. 

It  is  well  known  that  under  proper  condi tl  one  of  ph  asing 
between  mass  rate  and  ch.a  her  pressure,  spontaneous  oscillations 
in  the  florae  front  can  occur. 

The  procedure  adopted  was  to  assume  that  too  injection 
3yst«;  is  vibrating,  thus  establishing  synchronized  fuel,  flow 
modulation  and  then  to  derive,  on  tho  basis  of  a  simple  dynamic 
model,  a  formula  for  tiie  resulting  nrossurc  variation. 

Simplifying  assumptions  made  in  the  model  ire  oa  follows:- 

(a)  Shower-head  injectors  with  low  velocity  jct3,  so  that  under 
vibration  the  break-up  is  into  uniformly  sized  drops. 

One  thus  has  waves  of  irons  entering  tho  combustion  zone. 

In  practice,  of  course,  tho  jet  is  shattered  into  droplets 
and  discrete  drops  are  replaced  by  concentrations  (bunches) 
of  dronlots.  dynamically,  the  situations  arc  similar. 

(b)  The  cxidiscr,  say  liquid  oxygen,  is  vaporised  much  more 
repidly  ts.an  the  fuol.  Ponce  fuel  injection  only  is  considered 
of  significance  in  this  work. 

(c)  All  flow  is  axiaL. 

(d)  Only  the  longitudinal  mode  of  instability  will  be  specifically 
involved. 


The  full  Jc  tails  of  this  initi  >.L  .•|ori<  r.-ay  be  found  in 
rof  ,8.  “ha  eioiel  was  a  very  crude  analogy,  in  which  the  pressure 
variation  following  the  injection  of  one  *wavo*  of  fuel  droos  into  the 
combustion  zone  is  equated  to  the  time  response  of  a  linear  second 
order  system  following  an  imoulou-tyie  disturbance.  Tho  equation 
for  ores  sure  as  a  function  of  tirae  vna 

— ps  V  ^ 

whore  p  s  instnn baaoous  combustion  pressure 
^  «  damning  ratio 

^  «  time 

^  a  constant  (a  function  of  tho  specific  chamber, 

fuel/oxidant  ratio  etc.) 

shore  salons  for^lo^K  in  terms  of  system  parameters  were 
derived.  These  parameters  wore  tho  combustion  tirae  lap  and  residence 
time  of  tho  gases  in  the  chamber. 

The  intersection  of  the  pressure  variation  for  successive 
•waves’  Leads  to  the  high  frequency  combustion  oressure  variation 
as  illustrated  in  Figure  1.  /arinti on  of  frequency  with  injector- 
ores  sure  drop  is  s.iown  in  Fig.2. 


Conclusions 

This  initial  crude  analysis  ’was  sufficient  to  demonstrate 
that  if  fu.  1  flow  modulation  cm  bo  produced  in  the  frequency  range 
known  to  be  associated  with  the  ohonononon  of  high  frequency  com¬ 
bustion  instability,  tnen  a  large  amplitude  pressure  variation  will 
result. 


It  thus  became  essential  as  a  next  sten  to  establish  whether 
or  not  injector  a  Late  vibration  could  produce  tds  flow  mod  til  at  ion. 

»s  mentioned  previously,  most  of  the  work  on  t  iis  grant  has  involved 
the  investigation  of  fluid  flow  thro’  injectors  in  the  presence  of 
mechanical  vibration.  Thin  -111  bo  reported  in  Part  2. 


Part  1,  Appendix  1 

Rocket  Engine  Noise  and  Structure  Response 

It  is  well  known  that  the  most  important  source  of  vibration 
in  a  rocket  in  a  static  test  or  during  the  powered  subsonic  part  of 
flight,  is  tlie  sound  field  of  the  rocket  engine. 

The  noise  of  the  engine  has  its  origin  in  the  very  intense 
turbulence  created  in  tne  shear  layer  when  the  fluid  exchanges 
momentum  with  the  atmosphere  and  has  a  much  greater  effect  than 
direct  vibration  transmission  frcm  the  engine. 

The  turbulent  sound  producing  gas  extends  downstream  from 
the  nozzle,  and  two  distinct  regions  may  be  identified  in  this 
extended  source. 

(a)  the  "near  field*'  region  within  a  wave-length  or  so  of  the  nozzle, 
including  not  only  outwardly  propagating  waves,  but  also  reciprocating 
motions  and  pressure  fluctuations. 

’b)  the  "far  field'*  where  only  radiated  sound  (representing  energy 
abstracted  from  the  jet)  is  involved  and  where  the  pressure 
fluctuations  fall  off  as  l/distance. 

Work  carried  out  on  the  near-field  effects  over  the  vehicle 
surface  ^  indicate  that  they  are  small  (at  least  for  frequencies 
above  500  cps)  and  so  from  the  point  of  view  of  structure  response,  it 

is  only  the  far-field  which  is  of  interest. 

(2) 

Lighthill  hr.s  '  emonstrated  that  the  radiation  of  soundby 
turbulence  is  caused  by  a  volume  distribution  of  quadrupoles  -  lateral 
or  longitudinal.  The  noise  field  along  the  vehicle  surface  will  only 
involve  the  longitudinal  quadrupoles. 


(3) 

Measurements  of  the  sound  radiated  by  rocket  engines  ,show 

that  tn*;  ratio  of  sound  power  to  mechanical  power  is  a  constant, 

independent  of  the  exhaust  Mach  number,  the  value  being  in  the  region 
—2  8 

of  10  .  In  fact  the  V  relation  predicted  by  Lighthill  reaches  a 

saturation  point  somewhere  between  the  jet  and  rocket  engines. 

The  dimensionless  spectrum  function  for  rocket  engines, based 
on  experimental  measurements  is  shown  in  Fig.3.  It  is  to  be  noted 
that  the  spectrum  is  wide. 

For  a  St  rouhal  msnber  of  unity,  with  a  medium  engine  of  3  feet 

diameter,  140,000  lbs  thrust  and  7,000fps  exhaust  velocity,  the 

frequency  would  be  2,300  cps  and  spectrum  function  A(f)  *  1.1. 

The  spectrum  function  is  proportional  to  the  power  per  unit-  bandwidth 

<••4.  *  2 

and  corresponds  to  a  spectral  density  of  5  x  10  (lbs/sq.ft)"  c.p.s. 

The  mean  surface  sound  pressure  would  be  about  50  lbs/sq.ft. 

The  total  noise  pov^r  in  the  above  unit  would  be  in  the  region  of  five 
million  foot-pounds  per  second. 
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Effact  of  Missile  Motion 

There  arc  two  known  major  effects: 

(a)  the  radiated  power  increases  with  the  rocket  velocity.  Lighthill 
has  shown  that  the  differential  power  radiated  in  the  forward  direction 
by  longitudinal  quadrupoles  is  increased  by  the  factor, 


1/(1  -  M)  where  M  =  rocket  Mach  number. 

The  engine  noise  field  along  the  vehicle  surface  will  tend 
to  become  very  large  as  M  aporoaches  unity.  The  increase  is  limited, 
hovever,  as  it  becomes  increasingly  difficult  to  excite  vehicle  vibrations 
as  M  approaches  unity,  because  the  effective  sound  wave  length  is  reduced 
by  the  factor  (1-M) . 

In  fig.  4  a  graph  of  M  versus  spectrum  function 
amplification  factor,  B(M),  is  shown. 

(bj  the  directivity  pattern  turns  more  towards  the  direction  of  motion. 

Certainly  tnen,  the  vehicle  motion  leads  to  an  appreciable 
increase  in  the  noise  field  at  the  structure  surface. 

Response  of  the  Vehicle  Structure  to  the  Noise  Field 

As  long  as  the  rocket  Mach  number  is  below  unity,  the  sound 
generated  in  the  far  field  (this  may  be  centered  40  ft  or  so  downstream 
from  the  nozzle)  can  reach  the  vehicle  surface. 

Only  large  structural  response  is  of  interest  and  so  two 
particular  cases  need  be  considered: 

(1)  Resonance  -  which  occurs  when  the  frequency  of  the  exciting 
force  equals  a  natural  frequency  of  the  structure.  The  noise  spectrum 
being  a  wide  one,  all  vehicle  resonances  can  be  excited. 

(ii)  Modal  Equivalence  -  which  occurs  when  the  sound  wave-length 
(or  speed)  is  equal  to  the  wave-length  (or  speed)  of  free  flexural 
vibrations  in  the  structure  ^ 

The  sound  speeds  involved  range  over  a  relatively  narrow 
band  from  the  sound  velocity,  c,  to  zero,  for  rocket  Mach  number 
ranging  from  0  to  1,  and  so  modal  equivalences  will  be  rare. 

The  largest  vibration  response  will  occur  for  those 
situations  where  both  resonance  and  equivalence  coincide,  and  thus 
those  resonances  with  a  wave-length  equal  to  a  noise  field  wave-length 
are  of  maximum  interest.  Such  occurrences  may  be  called  Coincidences. 

Crandall ^  has  developed  a  step-by-step  approach  for 
determining  these  most  important  structural  responses.  He  introduced 
circumferential  modes  to  describe  the  angular  variation  of  both  the 
sound  pressure  and  the  induced  vibration. 


» 


-Ki¬ 
ll'  w  represents  the  flexural  (radial)  displacement  of  the 
structure  surface,  the  vibration  acceleration  is  w,  and 

Jo 

where  (j)(y)is  the  acceleration  spectral  density.  The  displacement, 
w,  is  represented  by  a  Fourier  series: 

O 


where  the  co-ordinate  system  (x,  is  as  shown  in  Fig.  5,  and  n  is 

the  circumferential  order  number. 

The  n  »  0  mode  corresponds  to  expansion  and  contraction  of 
the  cylinder  cross-section 

n  =  1  mode  corresponds  to  translation  of  the  cross-section 

n  2  modes  correspond  to  bending  or  corrigation  of  the 
cross-section  perimeter. 

In  circumferential  modes  the  sound  pressure,  ^>  ,  is  written  similarly 

-p  6:  ,<£,*■)=  Z,  -Py/x,!)  ore  y^<f> 

y\~0 

The  task  is  to  express  the  acceleration  spectral  density,  m 
in  terms  of  the  mean  square-pressure  spectral  density. 

It  turns  out  that  for  pressurized  cylinders  there  is  only 
one  circumferential  modeof  importance  for  low  frequencies  -  the 


as 


n  =  1  mode. 

For  higher  frequencies  (  ^1000  cp3  say)  all  modes  are  important 

and  behave  identically. 

It  also  appears  that  coincidences  can  only  occur  for  vehicle 
Mach  numbers  up  to  about  0.7. 

Calculations  on  hypothetical  rocket  vehicles  of  the 
140,000  lb  thrust  category,  show  that  coincidences  can  occur  during 
the  early  part  of  the  flight  (or  in  static  test)  and  that 

(a)  the  frequencies  fall  in  the  range  from  1,000  to  10,000  cps 

(b)  rms  accelerations  up  to  40g  are  possible. 

(c) n=0,  1  &  ^  2  modes  are  possible. 

Luperi  and  Tick  (6)  have  recently  reported  work  on 
vibration  and  unstable  combustion  in  the  UR  99  engine. 

Acoustic  pressure  oscillations  induced  resonance  in  the  engine 
mounting  structure,  resulted  in  build-up  to  vibration  accelerations 
as  high  as  400  g.  Even  with  isolators  installed,  values  of  150  g 
were  obtained. 


Associated  with  these  vibrations,  a  great  increase  in 
the  magnitude  of  the  combustion  prossurs  oscillations  i tM* 
This  wa3  ascribed  by  the  authors  as  possibly  being  connected 
with  the  effect  of  mechanical  vibration  of  the  injection  system. 


It  is  contended  in  this  work  that  the  initiating  mechanism 
for  all  forms  of  high  frequency  instability  in  liquid  fuel,  rocket 
engines,  can  be  traced  to  the  massive  mechanical  vibration  of  the 


injection  system  driven  by  the  intense  sound  field  arriving  at  the 


vehicle  surface. 


(fi'e-i*  C'v'V  ^ 
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u^S^4°Sls 

JL  •  the.  expansi  on  of  the  jet  surface 
|3  *  a  constant 
£  --liquid  density 

-  surface  tcnsi.n 

-  position  of  the  jet  surface 
«  wave  mmber  of  disturbance  on  jet 

-  wave-length  "  "  "  " 

»  exponential  build-up  factor 
»  liquid  viscosity 

=  acceleration  due  to  gravity  (32.2.  ft/sec  ) 
«  ordinate  along  jet  axis 

*  jet  raiiu3 
»  velocity  potential 
«  axial  jot  velocity 

-  radial  »  " 

«  constant 

»  modulation  velocity 

*  background  modulation  velocity 
«  angular  velocity  -  rrr-f 
«  frequency 

«.  potential  energy  of  jet  per  unit  length 

-  kinetic  energy  of  jet 
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Tha  authors  have  reported  initial  experimental 
work  on  the  effects  of  mechanical  vibration  on  the  break-up 
of  a  laminar  cylindrical  water  jet  in  air. 

Results  were  compared  with  the  1st  order  theory 

developed  by  Lord  Rayleigh  ^ .  He  had  shown  that  a  cylindrical 

liquid  jet  could  be  triggered  into  instability  by  small 

rotationally  symnetric  disturbances,  provided  the  wave-ler.gth 

on  the  jet  was  greated  than  the  jet  circumference. 

(3) 

C.  Weber  carried  out  a  similar  analysis,  and 
extended  it  to  include  the  effect  of  viscosity  and  an  estimate 
of  aerodynamic  effects  at  t^e  higher  jet  velocities. 

(above  about  15  metres/sec) . 

No  attempt  was  made,  however,  to  specify  the 
physical  nature  of  the  initiating  disturbance,  nor  to  extend 
the  analysis  to  the  case  of  finite  driving  amplitudes.  In 
this  work  both  of  these  points  will  be  covered.  The 
equation  of  free  motion  obtained  by  Rayleigh  was  as  follows:- 

_ d-i) 

A 3  JL  Cc  en  then 

,  .  -i-WT 

_  »  0\  V  /.  —  \ 


'l  ~  u-2) 

where  q  is  the  exponential  build-up  factor. 

Weber  extended  this  analysis  to  cover  viscous 
effects  and  the  amplifying  aerodynamic  suction  effect.  He 
obtained  the  following  eauation  for  q:  z.  ,  .  (iV  .  a  \ 

For  relatively  inviscid  liquids  such  as  water  and 
paraffin  the  viscous  correction  can  be  igiored  and  by 
considering  jet  velocities  below  1500  an/sec  the  aerodynamic 
effects  are  small  ^ .  Moreover,  as  will  be  seen  later,  in  the 
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presence  of  finite  initiating  amplitudes  Weber's  aerodynamic 
treatment  (only  applicable  to  sinusoidal  corrugations)  would 
have  to  be  replaced  by  a  more  appropriate  analysis. 

The  break-up  and  periodicity  of  liquid  jets  has 
been  studied  again  recently.  This  stimulus  is  duo  to  possible 
connection  with  the  phenomenon  of  combustion  instability  in 
liquid  rocket  engines.  It  was  found  ^  ^that  the  Rayleigh- 
V/eber  capillary  type  instability  could  be  driven  by  subjecting 
the  liquid  jet  to  sound  fields.  Very  small  pressure  variations 
and  low  frequencies  were  used  (lo,5  psi  and  200  cps)  and  the 
results  were  unspectacular  and  limited  to  low  velocity  laminar 
jets. 

Reba  and  Brosilow  ^  ^ ,  however,  extended  this  work  to  large 
pressure  variations.  The  jet  entered  a  chamber  which  effectively 
had  a  reciprocating  pi3ton  at  one  end.  The  pressure  downstream 
of  the  injector  was  varied  from  5  to  30 %  of  the  steady  pressure 
at  frequencies  up  to  3,000  cp3.  Pronounced  periodicity  was  induced 
on  the  jets  (oil  and  water) .  Considerable  quantitative 
investigation,  such  as  the  effects  of  oscillation  frequency 
and  amplitude,  injector  hole  size,  etc.  was  carried  out. 

Reba  and  Brosilow  postulate  a  mechanism  based  on 
flow  modulation  caused  by  the  pressure  variation  dovtistream 
of  the  injector  hole.  They  ascribe  the  corrugations  cn  the 
jet  as  due  .purely  to  the  resulting  "bunching"  action  of  the 
liquid  in  the  jet.  (bunching  is  a  term  used  to  denote  the 
radial  velocity  of  liquid  due  to  the  relative  velocity  of 
adjacent  particles  in  the  jet). 

Arbitrarily  they  ignored  capillary  action,  but  for 
the  low  initiating  amplitudes  obtained,  as  evidence  by  the 
sinusoidal  type  of  biild-up  on  the  jet  surface,  it  was  unjustified. 
This  is  characteristic  of  predominant  capillary  action. 
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Their  analysis  considered  ‘he  effect  of  the 
resistance  and  inertiance  of  the  fluid  in  the  orifice  cavity 
on  the  flow  modulation,  and  analytical  formulae  predicting 
jet  length  to  the  "first  thickening"  were  derived,  again, 
this  distance  criterion  is  hypothetical  as  'bunching’  is 
a  continuous  process  (till  do-bunching  starts). 

It  will  be  seen  later  that  Reba  &  Brosilow  would 
have  been  more  correct  in  their  assumption  of  predominant 
"bunching"  at  much  higher  velocity  modulation  values.  It 
is  doubtful,  nowever,  if  such  values  could  be  achieved  by 
reasonable  pressure  variations  downstream  of  the  injector. 

It  is  the  pressure  upstream  which  fixes  the  discharge  velocity 
and  the  downstream  pressure  waves  must  propagate  upstream  and  in 
doint  sc.  their  effect  is  onormously  attenuated. 

In  this  respect  injector  vibration  is  more  efficient 
in  that  it  appears  to  induce  pressure  variations  directly  upstream 
of  the  injector  hole  (see  Section  11). 

Thus,  the  only  theoretical  analysis  cf  significance 
in  this  field  has  been  that  of  Rayleigh  and  Weber  and  is 
restricted  to  1st  order  effects  and  infinitely  snail  initiating 
amplitudes .  With  the  latter  restriction,  of  course,  the  flow 
modulation  merely  serves  as  a  trigger,  liquid  bunching  is 
negligible,  and  caoillary  type  instability  predominant. 

It  was  the  objective  of  this  work  to  praiuce  a 
theory  which  (a)  will  cover  second  order  effects  (i.e.  will 
extend  to  appreciable  corrugation  sixes);  (b)  th am  will  include 
the  capillary  effect  and  the  'bunching'  effect.  The  latter  must 
be  considered  when  the  initiating  amolitudo  (velocity  modulation) 
is  of  finite  magnitude  and  will  became  predominant  at  large 
initi ating  amolitudes . 


The  first  analysis  (Section  111)  to  be  presented 
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will  merely  extend  the  Rayleigh  theory  to  Second  order  effects, 
but  will  iucluie  the  velocity  modulation  us  the  initiation  for 


the  instability  (that  is,  restriction  to  very  sman  velocity 
modulation  is  made). 


The  second  analysis  not  only  extends  the  treatment 
to  second  order  effects,  but  includes  the  'bunching'  effect 
and  so  attempts  to  bridge  the  transition il  gap  between 
preieminant  capinary  action  and  predominant  liquid  'bunching* 
action. 
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Section  2,  Experimental  Background 
# 

The  min  experimental  results  reported  hero  are  an 
extension  of  the  work  previously  reported  by  the  authors  ^ . 

Briefly  it  covers 

(a)  the  results  obtained  at  much  nigher  vibration 
acceleration  values  (up  to  above  300  g) . 

(b)  the  rope'^ng  of  break-up  length  versus  initiating 
vibration  amplitude  measurements  for  paraffin  and  coma  orison  with 
the  water  results.  It  is  assumed  that  the  induced  velocity 
modulation  is  proportional  to  vibration  amplitude. 

(c)  the  detection  of  pressure  fluctuations  in  the 
cylindrical  chamber  behind  the  inject  r. 

details  of  m„st  of  the  aoparatus  us~d  has  been  giv<.n 
in  Reference  1. 

Previous  results  my  be  summarized  as  follows 

For  a  lamin-a*  water  jet  subjected  to  modulation 
driving  amplitudes  corresponding  to  vibration  acceleration  values  of 
up  to  20  '  g) ,  the  graph  of  breakup  length  (or  time)  fallowed 
reasonably  closely  that  predicted  by  the  Rayleigh- ^bb~r  tneory. 
However, 

(i)  using  an  experimentally  determined  discharge  coefficient  to 
get  the  true  initial  jet  radius  following  contraction,  estimates 
of  the  buil  i-up  factor,  q,  indicated  a  value  of  ab -ut  30C.  This 
is  significantly  higher  t.nan  the  value  of  225  expected  on  the 
Rayleigh-Weber  theory. 

(ii)  the  graphs  of  instability  build-up  versus  jet  length 
(or  time)  rev-.  il=d  c'nsiderable  departure  frar.  the  purely 
exponential  characteristic  inherent  in  tne  Rayleigh  analysis. 

As  mentioned  before,  the  slightly  peaked  appearance 
of  the  corrugations  at  vibration  acceleration  values  of  about  20  1 g’ ^ 
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was  interpreted  as  an  indication  that  velocity  modulation  is 
at  work.  At  low  accclorati. in  value*  this  mo  halation  is  very  small 
and  acts  only  as  a  trigger  for  the  Rayleigh  capillary  type 
instability.  *t  higher  valuo3  and  thus  stronger  modulation,  liquid 
bunching  along  the  jet  would  become  appreciable. 

To  investigate  this  further,  vibration  values  up  to 
about  250  ' g*  were  aopiied  to  the  cavity  behind  the  injector. 

Typical  results  for  water  and  paraffin  are  shewn  in  Plates  11.  and 
111.  The  jet  velocity  was  45  fps,  and  tne  vibration  frequency  2,000 
cps. 

Disc  formation  cluse  to  the  injector  is  the 
extraordinary  result.  The  paraffin  discs  are  more  affected  by 
aerodynamic  effects,  being  bent  and  atomized  round  the  edges. 

This  is  probably  lue  to  paraffin  h  iving  a  much  smaller  surface 
tension  than  water. 

The  bulges  between  the  discs  (in  the  case  of  water) 
is  due  to  the  build-up  of  harmonics  of  the  aopiied  vibration  on  the 
cylinders  joining  the  discs. 

The  disc  formation  is  striking  confirmation  of  the 
valocity  modulation  mechanism. 

It  is  obvious  that  the  Rayleigh-Weber  theory  could  not 
b-  expected  to  cover  this  situation. 

k  satisfactory  theory  for  finite  initiating  amplitudes 
(velocity  modulation)  must  include  both  capillary  and  bunching 
effects . 

Graphs  of  break-up  length  versus  log  (vibration 
amplitude)  f  or  paraffin  and  water  j-ts  at  a  velocity  of  1.5 
metres/sec  are  shown  in  Fig.l.  d'-tc:  the  ordinate  scale  is  in 
accole-rcmeter  volts.  This  is  prop-rti onol  to  vibration  amplitude, 
in  turn  pr^po^ti 'nil  to  acceleration  value.  It  is  reasonable  to 
assume  (see  discussion  lat<_r)  that  the  velocity  modulation 
amplitude  is  linearly  related  to  the  1 g*  value.  Experimental 
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points  weru  token  up  to  -about  50  '  g’ s  of  vibration  (  about  3 
volts  accelcraiictor  output) . 

Thu  curvos  for  water  md  paraffin  are  Vc. ry  similar. 

The  cross-over  point  W'uli  inlicatc  t.o.ot  for  a  given  vibration, 
a  larger  velocity  modulation  value  is  induced  in  paraffin.  As 
the  effect  is  probably  closely  connected  with  the  liquid  elastic 
properties  at  high  frequencies,  one  might  expect  this,  03 
paraffin  has  a  lower  relaxati  in  frequency  than  water  and  so  a 
higher  shear  modulus  of  -.lasticity. 

Thu  curves  theoretical? y  should  approach  infinite 
break-up  lengths  at  zero  amplitude,  but  very  small  randan  effects 
(pressure  fluctuations  in  the  injector  cavity)  are  alv/ays  present 
and  cause  the  finite  limits  observe  l  in  practice. 

It  is  fairly  obvious  that  there  is  no  unique  exponential 
buili-up  factor  1  g*  .  However,  by  restriction  to  the  range  below 
about  1.0  volts  (or  20  ’  g1 )  an  approximate  q  factor  can  be 
extracted  and  this  is  in  effect  what  was  d;ne  in  cur  previous  work. 

nn  attempt  was  ar-’e  t,  isolate  the  bunching  effect 
(dependent  on  frequency)  from  the  c-piliary  effect  ( do  indent  ...n 
wave-length) .  .a  series  of  break-up  length  versus  vibration 
frequency  measurements ‘was  carriu  !  out,  keepinr  the  wave-length 
on  the  jet  constant  by  variation  of  jet  velocity.  The  resultant 
graph  (for  water)  is  shown  in  Fig. 2  (break-up  time  has  been  used 
here).  It  is  s-  on  that  ob'-ve  ab  ut  l,0eC  cps  the  break-up  time 
decreases  almost  linearly  with  frequency  (cr  acceleration  value,  as 
the  amplitude  was  also  hull  constant).  3elcw  about  1,000  cps 
vibration  wouli  have  little  effect  on  break-up.  It  is  o-.-ssible  that 
this  corresoonds  to  a  first  relaxati  -n  frequency  for  water. 

3y  the  use  of  very  small  velocity  m^iulition  values  it 
should  be  oossible  to  simulate  the  infinitely  small  disturbance 


1 


* 


* 


-  22  - 

region.  This  was  lone  by  using  a  frequency  of  1,000  cps,  a  jot 
velocity  of  100  cm/sec  and  ' g'  from  4  to  20.  The  resulting 
curve  is  shown  in  Kir. 3,  It  is  3een  that  a  good  straight  lino 
is  obtained.  A  q  value  of  293  was  obtained  from  this  graph.  It 
will  be  seen  later  (Section  3)  tnit  the  difference  between  this  and 
the  value  expected  using  Rayleigh's  analysis  is  explained  by  baking 
2nd  order  effects  into  account. 

High  ,/e locitv  Jets 

Alti  ugh  the  theoretical  analysis  will  be  restricted  to 

laminar  jets,  one  would  exooct  that  tne  strong  liquid  bunching 

produced  with  very  high  1 g*  vibration,  should  have  an  effect  on 

high  velocity  turbulent  jets.  This  of  course  is  the  case  most 

relevant  to  rocket  engine  combustors.  That  this  is  true  can  be 

clearly  seen  from  Plates  IV.  and  V.  These  photos  are  of  water 

jets,  from  a  0,25  an  orifice.  The  velocities  were  90  and  145  fps 

respectively,  with  vibratirn  acceleration  of  about  200  '  g'  and  at 

a  frequency  of  4-500)  cos.  The  iisc  f  nnition  is  still  very  evident, 
1^800) 

but  ligament  formation  with  subsequent  break-up  into  droplets, 
tends  to  cloud  the  iiscs  from  view. 

To  prove  that  mass  flow  modulation  was  occurring  even 
with  high  velocity  jets,  a  focussed  light  beam  was  passed  tnrough 
the  jet  and  on  to  a  photo-electric  cell.  The  output  from  the 
ohoto-cell  was  fed  t:  a  T'arconi  wave-form  analyser  uai  the  spectrum 
in  the  vicinity  of  tiie  aopliud  vibration  frequency  examined  carefully. 

Up  to  the  highest  pressures  possible  (about  150  p.s.i. 
yielding  a  jot  velocity  of  roughly  140  fos),  a  sharp  peak  in  the 
analyser  outout  was  found  at  the  apolied  frequency.  Th^  results 
are  illustrated  in  Fig. 9  (background  noise  has  been  subtracted  cut 
and  was  always  at  least  an  order  of  magnitude  below  tbv.  peak) , 

One  may  conclude  that  vibration  will  induce  appreciable  miss  flow 
modui  ition  even  with  high  velocity  turbulent  jets. 
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Injector  Cavi ty  Pressure.  Fluctu  \ti  ;ns 

Rcb'i  mi  3r.r>iiow  revortei  tn  .t  they  detected  no  pressure 
fluctuations  in  the  cavity  behind  the  injector,  as  the  mechanism 
of  velcity  me  iul  .ti  ..n  w  uli  involve  pressure  fluctuations  upstream 
of  the  injector  Inlc,  an  attempt  was  made  to  detect  these  using 
a  high  frequency  Piece  transducer.  This  had  a  resolution  of  0.015 
psi,  an  unper  frequency  limit  ef  over  10,000  cps  but  was  designed 
for  a  pressure  range  <1  up  t.»  200  psi  (giving  a  3  volt  output). 

The  greatest  iifficulty  lay  in  the  fact  that  the 
transducer  was  sliehtly  sensitive  the  vibration  and  all  readings 
had  t .  be  subtracted  from  a  backgr  und  level  taken  with  no  liquid 
in  the  cavity.  hLs  m  iXimum.  vibration  acceleration  values  of  only 
up  t)  about  25  1  g'  were  safe. 

The  results  were  thus  net  conclusive  but  did  indicate 
that  an  effect  (or  signal)  was  obtained  only  when  tne  transducer  was 
mount-  d  on  th„  fr.-nt  plate  cf  the  injector  cavity  -  perpendicular 
tv  tne  direction  :f  vibrati-.n.  A  signal  to  noise  ratio  of  2:1 
was  obtained,  the  f  llowinc  numerical  results  being  recorded: 

Vibration  ’ g*  :  25 

"  frequency:  4,500  cps 
Output  si  mol  (fr-.i't  PieZv  amplifier)  =  0.032 
"  !l  frequency  =  4,500  cps 
The  signal  represented  a  pressure  variation  of  I  2  p.s.i.,  or 
-  %  of  the  static  pressure  in  thfc  cavity. 

Pr-  jw-ction  .,f  the  transducer  f  ace  1/4"  beyond  the  plate 
surface  resulted  in  zero  signal  again. 

The  production  of  shear  waves  in  the  boundary  layer  along 

a  plate  oscillating  in  its  own  plane  is  1  known  a. id  has  beer. 

(7'  (3) 

studies  by  Glauert  1 ,  Schlichting  ani  others. 

Further  work,  tne -re ti cal  and  experimental,  on  the 
effect  of  vibration  on  the  hydrodynamics  of  the  injector  ca/ity  will 
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be  well  worthwhile.  It  will  be  tentatively  postulated  here  that 
shear  waves  produced  by  the  vibration,  at  the  injector  plate,  can 
propagate  across  small  injector  noles  and  thus  cause  velocity 
modulation  in  the  emerging  jet. 

Section  3.  Theoretical  Analysis  of  Jet  If  stability  for  V^rv  Small 
Velocity  Modulation  Amplitudes 

Consider  a  frame  of  reference  moving  with  the  mean 
jet  velocity.  With  axial  symmetry,  the  two  co-ordinates  of  the 
system  are  r,  the  radius,  and  z  along  the  jet  axis. 

The  surface  of  the  jet  at  any  instant  t,  is  given  by 

r  -  ^  *  oL  - (3-D 

where  a  is  the  radius  of  the  undisturbed  jet,  ani  i=  ^(z»  0  is 
the  build-up.  It  will  be  considered  initially  that  is  small 
with  respect  to  a. 

As  the  motion  is  irrotational, 

^<P~  o 

w  .ere  ^  is  the  velocity  potential  and 
u  *  axial  velocity  = 
v  =  radial  velocity  = 

'br 

The  boundary  conditions  which  <f>  must  satisfy  are: 

(i)  the  free  surface  condition,  i.e.  a  particle  on  the  surface 
will  remain  there. 

Consider  a  quantity, 


S  »  a  +Ji  -r.  This  is  a  quantity  which  varies  from 
particle  to  particle,  but  for  a  particle  on  the  surface  £  =0. 


Thus,  on  the  surface, 


0,  which  is  differentiation 


following  the  motion  of  the  fluia.  It  can  be  written  in  partial 
derivatives  as, 

S  .  IS  j.  IS  -At  . 

^  4_  vj‘  3E  +  vr  -3GE  =0 

or 

IS  'iS  '.S 

■bt  +  “  +  =° 


vi 


\  -r 

u  • 
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In  terms  of  i  this  yields  the  equation  corresponding  to  the  free 
surf  ice  boundary  condition. 

V  ,  ^ 

-4-  ^  -  0  - (3.2) 

(ii)  Bcrnouilli' s  equation 

If  p  is  tnc  pressure  in  the  jet,  and  ^  is  the  liquid 
density,  t.-.ir,  equation  is 


^  +  % 


-  (3.3) 


On  the  free  surface,  a  balance  of  forces  gives 

\  =  £^!i  +T/efe 


(3.4) 


where  R  is  the  jet  radius  of  curvature 
T  "  "  :I  surface  tension 

On  the  free  surface,  therefore,  the  following  equation  holds: 


V 


4" 


-(^)+(l4)]  +•  t4r 


constant 


-  (3.5) 


This  is  the  second  boundary  condition  equation.  It  can  be 
expressed  in  terras  of ^  by  means  of  a  tedious  calculation,  giving 


1 


-  -A .  - i. 

~  (  fe,  /  -r  _"C 


Vi 


^  ^  * 


+r  + 


Terms  ini.  have  been  neglected. 
First  ^ooroximation 


x? 


V 


1 3.6) 


.Neglecting  terms  ini  ,  we  assume  that  ^  has  the  fora 

=  fl.gV'hXj.kg _ 


and  ti  is  loads  tc  the  original  Rayleigh  expression  for  the 
exponential  build-up  factor  q. 
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Second  Approximation 

In  this  cise,  powers  of  &  must  not  be  neglected. 

_ t 

To  extend  the  analysis  to  second  order  terms,  (p  must  include 
a  term  corresponding  to  the  residual  term,  — *  — —  ,  in 

H  ^ 

the  boundary  conditions.  So  for  a  second  order  analysis  we  assume 

<*  =  $,+  &  =  f\ex*zikvTJiS*r) 

-  (3.8) 

JL  *c  *in^  ^  ts  1  constant 4 

.'ow  r  a  a  ,  anl  writin.'  the  boundary  conditions  to  a  second 
order  lpproxim iti or. ,  aid  taxing  t'.e  special  case  of  k  i  »  0,5 
(centre  of  the  Rayleip.h  spectrum) ,  lea  Is  to  the  following  expression 
for  (£\  , 

eiWUi)_ 


(3,9) 


where  q  -  Rayleigh  build-un  f  actor 

*,?:«  *VG-iv) .  £• 'r 


o-s* 

The  equation  I  or  j[^  tnen  becomes 


'V 


(3.10) 


A  value  for  the  constant  may  b-.  obtained  by  applying  toe  boundary 
condition  (or  initial  condition)  at  trie  orifice.  T;. is  will  be-  a 

I 

small  velocity  modulation  It  ,  and 


I 

5=  modulation  velocity  =  ~ 


- (3.11) 
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It  is  issumed  that.41,  is  known. 


Hence 


i  'x>^  4^*:)  Yve*“~  Y»J 

\  -  (3.12) 

A35=U.i  one  can  then  either  compute  the  build-up,  , 
versus  time,  ~r  the  profile  on  the  jet  it  successive  intervals  of 


time . 


liquation  (3.21)  can  be  reduced  to  the  following 


anoroximate  frrm: 


-  (3.13) 


A  typical  result  for  the  profile  on  a  water  jet  from  a  0.12  cm 
injector  is  si.own  in  Fig. 4.  Tho  jet  velocity  was  1000  cm/sec 

and  it  ms  assumed  t:iat  the  applied  vibration  was  producing  a 

v 

velocity  modulation  of  At  =  5  cm/sec.  Fig. 5  shows  values  of 
J ^  versus •<^'U(t  measured  from  t  =  C  at  the  orifice).  As  the 

build-up  will  be  still  closely  exponential  we  can  write 

/r  ^-6*}  ‘ic\ 

js  constant -Cd. 

•where  )  is  the  independent  variable. 

From  the  graph  a  value  of  K  =  1.5  was  determined  and  so  an 
effective  * qT  factor  of  330  was  arrived  at. 

This  value  agrees  quite  -well  with  the  experimental  estimates  cf 

(I) 

q  reported  in  the  previous  paper 


It  is  interesting  to  note  that 

\jV/2t«  =  V  T/2C-a 


and  so  has  the  dimensions 


of  velocity.  For  the  injector  used  and  with  water  as  the  liquid, 
this  characteristic  velocity  turns  out  to  be  about  50  cm/sec. 


Z  JU 


v. 


T 


8 

\TT 


(3.14) 
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The  series  represented  by  equation  (3.20)  becomes  divergent  for 

A 

an  value  over  one.  So  the  analysis  woula  break  down 

in  this  case  for  a  value  of  velocity  modulation  above  about 
12  cm/sec.  This  is  all  rignt  within  the  initially  assigned  limits 
of  very  small  initiating  modulation  amplitudes.  At  higher 
amplitudes  the  analysis  would  have  to  be  replaced,  as  liquid 
bunching  as  well  as  capillary  instability  must  b^-  taken  into 
account. 

It  is  felt  that  the  above  analysis  is  an  improvement 
on  the  original  Rayleigh  analysis  as 

(a)  the  physical  nature  of  the  initiating  disturbance  is  identified 
and  incorporated  in  the  theory. 

(b)  the  analysis  has  been  extended  to  the  second  order  and-  thus 
can  be  expected  to  be  accurate  for  appreciable  sixed  corruptions 
on  the  jet. 


! 

\ 

I 
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Section  A.  A  Theory  for  Jet  Break-up  in  the  Case  of  Large 
Velocity  Modulation  Amplitudes. 

The  velocity  moltilatien  in  this  cose  not  only  triggers 
the  capillary  instability  but  causes  liquid  bundling  along  the 
jet  (a  forced  radial  compounded  velocity  is  introduced). 

Bunching  Theory 

As  a  preliminary  step  the  bunching  mechanism  will  be 
considered  and  mathematically  defined.  The  jet  car.  be  considered 
as  a  system  for  particles  without  surface  tension  or  interaction. 

The  velocity  of  partied es  emerging  from  the  orifice 

at  time  t  =  t  ,  is  written  as 

o' 

AV  ,  ~  JXC  -ir  wtA  Swv  JLaj  - - (4.1) 

whore  t  is  the  time  original  at  the  nozzle. 

jU  o  is  the  steady  velocity 
AX  is  the  modulation  amplitude. 

The  time  of  arrival  of  a  particle  at  a  plane  z  is  given  by 


'O 


±,  =r  ± 

*  4*  1 *4~  i  — *ic~ 

•  •  A  —7  -'o  \  rV* o 


V 


U.  0  ("  (4-  ix  ^  ^ 

~  ,  -XAo  ^ 

As  the° flow  is  continuous,  if  nQ  is  the  particle 


(4.2) 

(4.3) 


current  at  the  orifice,  then 


-eUc  t 

o 


-Xxi  \  -U  ^ 

- *Z - *  CcJl  ,kKlX0 


-  (A. A) 


V\c 


1 

•  • 


T\ 


»  I  — - -  (  — o  )<^r$AJJ±0  u*5) 

If  the  particle  density  is  constant,  then 


Hi 

TV 


X 


/ 


T, 


where  }  T~,  ,  are  the  radii  of  the  jet.  The 
following  expression  fer  the  radius  of  the  jet  at  -any  time  t  (in 
the  orifice  frame  of  reference)  is  obtained: 
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nr,  =  l  4 


jlA 


&rs  vj  tc 


(4.6) 


It  may  be  sometimes  more  useful  to  consider  the  various 
harmonics  sen-,  rate  lj»,  is  the  above  expression  assumes  that  they 
are  equally  attenuated.  So  we  expand  outY^,  of  equation  (4.5) 
in  a  Fourier  series  in  the  angle  #  Considering 

the  first  harmonic  only,  leads  to  the  following  exoression  for  n: 

w,=  w.[  I  +  xl, 

-  va«.f  \+\Z(  - u.7) 

f,  -  m 

where  ^  —  _u0dt  -=  4TT/  ^  ,  and  the  time  scale  hos  boon 


,  and  the  time  scale  hos  boon 


change"!  to  that  of  the  moving  jet. 


In  Fig.  6  the  build-up  profile  for  the  case  of  all 
harmonics  present  is  shown.  It  would  indicate  that  with  pure 
bunching  (non-interacting  particles),  discs  would  be  formed  along 
the  jet.  As  will  be  seen  later,  this  prediction  was  strikingly 
confirmed  for  the  case  of  very  high  *  g'  vibraticn.  Under  these 
conditions  appreciable  velocity  modulation  is  induced  and  so  trie 
‘bunching*  action  predominates  over  surface  tension  effects. 

The  non-interaction  condition  on  which  the  above  equation  (4.7) 
is  base!,  -./ould  then  naturally  be  a  more  realistic,  approximation . 
Jet  Anal,  sis 

The  Hayleigh  procedure  was  followed  closely,  but  the 
surface  position  equation  was  modified  to  include  a  ’bunching7 


term: 


T  -=^G  4  ^  fo) 

■AA  a 


- (4.8) 


This  lecads  to  changes  in  both  the  kinetic  and  potential  energy 


terms.  The  kinetic  energy  must  now  include  the  transverse 


velocity  component. 
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The  expression  for  a  now  becomed, 

l./s'1  .  X 

i  _L  / a  ->  -A  4o 

<tx.  —  ,-<t s  +4  a  H - , — ~ 

-u.0 


_/*  l.  cA 


and  so 


^  J.  i  2- 

'^•o  —  —  4- 


I  -f- 


-1-  ,t  2-t 
ft  _M  j^J  \ 

U ■  Xk^  J 


-(4.10) 


■iotv.:  the  cosine  term  in  the  ’bunching1  expression  is  taken 
as  unity,  so  t.hnt  .nty  fha  px.'ofiV-  ->f  movun  v  m  t,he  build— up 
will  be  obtained. 

The  potential  energy,  P,  per  unit  length  now  is 
given  by  the  following,  relation: 


T  =  (i  %?’- 0  +  iAV 


+ 


M  MX 


(i+axY 


- (4. 11} 


!4~ 


0  2  ,**-  T-> 

U-  XX  c  J 


Kinetic  Energy  Term 

If  0  is  the  velocity  potential  then  the  axial  and 
transverse  velocities  are  ;dven  by 

jU  = 


Now 

where 


~ - (4.12) 

- (4.13) 


p.  *.<.  .  y  .. 

K,  the  kinetic  energy  .x.r  unit  length,  is  given  by 

/•  TtT-n.  r  <kr 

K=i  ^  (jJUr-f4-foU$r))]rMz„(i.u) 

-%-0  Xo  -  ®  A  7,  , 

and  using  the  relations  (4.12)  and  -p  :=o  xx  t  w<-  £vf/ 

K-IV  Aft  €  ^  [(%(  iJk'v))  +(l0(j-  A^)Jti4£(4.i5) 

Xo  ® 

The  Bessel  functions  are  now  expanded  out  as  series, 
terminating  those  at  the  third  and  second  terms,  respectively. 

r .  i+w 


The  final  expression  for  X  is,  «  ^ 

rrirVjfeftV^Fll 


<k 


4 


4- 


C  -2 

w  T— 


+  *- 


DENSITY  (relative) 
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Using  the  Lagrange  equation  for  a  conservative  system, 


IK  .  *? 


vk 

and  considering  maxima  only  ( ^9  —  |  \  •  the  following  equation  of 


VJL.  +  'U 


O 


I 


(> 


,1  't 

>u 


■Jb') 


4-jO 


4! 


T.'4  -M 


.17) 


It  must  bu  acted  that  in  ioriving  this  equation  all 
the  harmonics  wer-;  assumed  to  be  present  and  so  it  will  be 
strictly  only  valid  for  the  conditions  of  strong  bundling  (very 


high  * g*  high  frequency  vibration),  .-or  the  transition  region 
where  capillary  and  buncoing  action  are  roughly  of  the  same 
effectiveness,  an  analysis  based  an  equation  (4.7)  would  probably 
be  more  aporepriate. 

Method  of  3q Luti on 

Equation  (4.17)  must  be  solved  for  as  a  function  of 
time.  I  his  can  only  be  done  .numerically.  .\'cw  *V~  is  a  function 
of  both  ^  and  7^7  and  so  one  is  forced  to  compute  for  a  small 

interval  of  time,  comoute  T”  and  then  use  this  value  in  the  corn- 

put  ati  m  of  4^  at  the  end  of  tho  next  time  interval. 

Tijc  eth  >ricr  Runge-.autta  technique  for  inte, grating 
differential  equati ms  was  found  to  be  the  most  suitable  for  this 
pur^  se. 

v7hu n^,  is  iotc  mined  for  a  step,  the  new  y*  value  is 
computed  and  this  is  usu ,  in  equation  (4.S)  in  computing  ^  at  the 
end  of  the  f  ull 'wing  step. 

This  proce Lure  appears  to  give  satisfactory  results. 

The  programme  was  written  in  Fortran  and  executed  cn  t.ne  IBM  1620 
digital  computer.  Bessel  functions  were  evaluated  in  the  programme 
using  the  continued  fraction  approximation.  This  increases 
computing  time  greaoly,  however. 


Th.re  is  nc  U.u->reticol  analysis  available  for 
predicting  the  velocity  halation  induced  by  a  particular 
vibration  l  r\  qu  nev  and  1  o'  v>lue. 

The  ra®  ef  v  lues  was  chosen  ,Nn  the  basis  of 
the  deserved  pr^ssur-.  vari'tions  in  the  cavity  as  re aorted  in 
Section  2. 

l 

For  a  IXC  an/ sec  jot,  Xi  values  from  ab'ut  3  cm/sec 
to  30  cm/sec  were  used. 

This  w.uld  correspond  r'mgnly  to  acceleration  values 
from  about  5  t^  50  ' g* .  Modulation  values  less  than  about  5  cm/sec 
woull  bv  only  caa able  of  initiating  capillary  type  instability. 

Ih'_  re  -d-.n  front  5  to  X  cm/sec  is  transitional,  where 
capillary  action  and  bunching  are  almost  equally  effective. 

In  the  region  over  30  cm/ sec  liquid  'bunching*  becomes 

predominant. 

The  break- uo  length  criterion  has  been  interpreted  as 
that  mint  cn  the  jet  wo --re  the  maxima,  or  ,  has  reached  a  value 
equal  t-  the  original  jet  radius.  It  has  significance  at  lew 
or  intermediate  modulati  n  values,  where  bunching  \ction  is  small. 
However  with  strong  bunching  it  loses  its  significance  as  a 
criterion.  Disc  formation  corresponding  to  ^  values  of  many 
times  the  original  radius  is  very  rapid,  but  actual  break-up 
is  delayed  till  capillary  action  thins  out  the  liquid  cylinder 
joining  the  discs. 

Thus,  caparisons  of  experiment  and  theory  on  the 
basis  of  break-up  measurements,  is  only  oossible  in  the  low  and 
intermediate  'bunching*  regions. 

In  practice,  moreover,  small  random  fluctuations  in 
pressure  are  always  present,  and  so  one  cannot  have  zero  velocity 
modulation  as  in  theory.  To  make  the  above  theory  more  realistic, 
then,  a  small  constant  modulation  term  ,  was  included  in  the 
formula,  whore  lx  *=  J1  v  i-  c.  0 


In  the  1000  an/ see  c-ase  taken  as  an  example, 
w-'i3  given  .  v  due  of  2  ar/'sec. 

Comparison  of  Theoretical  and  Experiment'd  Results 

The-  ab-ve  analysis  was  appli-i  tz  the  case  of 
Jet  -at  1000  cn/sjc.  v docity,  with  an  applied  vibration  f 


6 


o 


a  water 
ruquency 


of  2000  cos.  The  wav. — length  on  the.  J-t  was  thus  0.5  ems 
(n^ar  th:  mi-1  11..  f  tnc-  Haylciah  spectrum). 


Oraphs  f  buiLi-uo  .profile  of  maxima)  -n  tr«e  .jet 
versus  time  are  sr.-own  in  tig.  7  for  2  cases.  Tht,  curve  for 
low  driving  amolituies  (weak  velocity  r.  iul  iti  n)  indicat-s  a 
fairly  rapid  Initial  build-up,  then  a  slow  increase  uad  final 
rapid  ’break-up'. 

The  experimental  result,  using  the  travelling 
microscope  technique,  is  shown  to  resemble  it  in  outline  (of  course 
iiie  initial  region  cannot  be  checked  as  the  .-.axima  are  too  small 
to  be  observed) . 

The  curve  for  higher  velocity  modulation  (30  cm/ sec) 
shown  the  raoid  build-up  waich  results.  The  break-up  criterion 
corresponds  to  an  orlinatv,  value  rf  .06,  and  has  no  significance 
in  this  case  as  mentioned  .above.  Plate  11.  shows  the  result, 
observed  for  strong  bunching.  This  was  a  water  jet  subjected  to  a 
vibration  acceleration  v alu<~  of  ab-'ut  60  '  g'  at  3,000  cps.  The 
rapid  disc  fomati'n  is  characteristic  of  'bunching'  in  which 
all  the  hanr.onics  'ire  present.  The  disc  reached  a  certain  size 
and  final  br-nk-u?  occurs  furtaer  down  the  jet  and  is  due  to 
capillary  action  on  the  cylinder  connecting  the  Uses.  The 
bulge  midway  between  the  discs  indicates  the  build-up  of  other 


harmonics . 


The  observation  of  disc  formation  is  striking 
confirmation  tr.at  the  driving  force  is  velocity  modulation. 


THEORETICAL  PREDICTION  OF  AMtJTUDE  BUILD-UP  VERSUS 

V  ■  JET  V£LOCITY  -  1,000  cu.Utc. 


FJG* 

WEAK -UP  LENG1 TH  Yk  VELOCITY  MODULATION  -  thcoaltical  - 


ACC£L£AOkt£T£R  VOlTS  (viamATiON  amplituoc) 


Plato  111.  sh  ws  similar  ri'nen  n.non  for  .  paraffin  jet. 


Most  of  til  experimental  me  asuronents  dado  on  the  jet 


•  1/ 1  ^  '•+  »■*  <*/»**«■»»»<•*  iri  Vav*  <%  f  i  ->n 

4V“UW  J.V  »  t  «  w»*  v  VI  b>\«u  V  J.  V*  .  ••• 


1  1  _ _ _  v  ** 

•O.XJ  U  i  i£-'  .i  V. 


thv-S-. ,  break-up  lengths  (based  on  the  radius  criterion)  for  a 
range  of  1  w  and  interne  ii  ate  moiul’.tion  values,  weru  computed 
fer  the  lO'JO  cm/soc.  j-.  t.  Vibration  amplitude  and  v^l^city 
modulation  are  considered  as  Proportional  parameters. 

Results  snown  in  Fig.  8  curves  (a)  and  (b) 
are  for  water  and  paraffin  jots. 

Tdu  th_  ary  predicts  little  difference  in  their 
behavieus  in  the  finite  modulation  velocity  region.  Curves  (c)  and 
(d)  are  exaerimcntal  ones  for  paraffin  and  water  jets  at  about  1200 
cm/ sec . 

The  experimental  curves  also  show  that  in  the  finite 
modulation  region,  the  behaviour  r  f  paraffin  and  water  is  indeed 
similar  (and  prosu.ably  this  will  be  even  more  so  at  higher 
modulation  valued). 

There  is  e.n.  ugh  agreement  between  tne  above  analysis 
and  experimental  results  U-  confirm  that  it  is  largely  correct. 

The  photographs  of  jet  break-up  -btained  by  Reba  and 
Brosiiow  and  Mi^sse  indicate  that  the  velocity  modulation  obtained 
by  then  falls  in  the  very  low  category,  covered  by  the  analysis 
in  Section  111.  Thus,  the  variation  of  pressure  i 7wnstre-am  of 
the  inject >r  is  not  an  efficient  method  of  oroducing  velocity 


modulation. 
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Section  Y.  Conclusions 

It  sterns  f  iirly  conclusive  tmt  the  rot  >.ti  -nally 
symmetric  disturbance  require  !  t  trigger  the  RoyL-igr.-Yubcr 
capillary  instability  on  liqui  i  jets,  may  lx-  ohysicolly  attributed 
t  ■  velocity  modulation.  This  w  uli  b-  induce  i  by  pressure 
vari  .tion  upstream  >.f  tiio  inject  r  thro 'll.  bn'U.r  rr  rmol 
circumstanc -s,  minute  ran  i  r,  oressuro  fluctu  .tions  will  always 
be  ore sent  in  the  cavity  behin  i  the  injector. 

The-  opplicati  n  if  .  :ecr  onical-  vibrati  in  in  the 
nppropri  ite  frequency  rmge  on  1  at  small  vibration  acceleration 
values  con  induce  such  minute  pressure  fluctuotions  and  trigger 
the  co’iilliry  instobility.  The  theory  exj  undod  in  .lection  3  is 
considered  ti  cover  tnis  cose. 

’nth  nigh  vibration  acceleration,  finite  velocity 
r.oduloti,  a  values  can  be  >r  - luc  ;d .  This  nr‘’iuces  on  opprecioblo 
'bunching*  effect  in  th<.  liquid  jet,  -end  if  strong  on  -u?h  it  will 
predoninote  -ver  the  c  * pillory  effect.  The  onolvsis  in  Section  k 
is  given  as  o  semi-quantitative  theory  of  the  phe lonenen  os 
observed  in  the  tronsiti >n  regi  n  between  the  nredontinant 
capillary  region  (infinitely  small  ' bunching' )  ond  the  oredominant 
'bunching'  region. 

iiowve r ,  o  crii.ori  n  tr  r-pl-ce  the  break-up  length 
-ne  (used  with  success  in  tiie  very  low  'bunching'  region)  is 
necessary  t  >  allow  satisfactory  c  n^aris on  bu tween  theory  and 
experiment . 

The  analysis  ire  iicts  that  under  strong  'bunching1 
conditions  (induced  by  high  frequency,  high  '  g'  vibration),  the 
behaviour  cf  various  liquids  will  be  very  similar,  Zxoeriment 
confirms  this.  Seme  variations  will  be  expected  due  to  the 
variation  in  c-lastic  oroperties  of  liquids.  For  a  given  vibration, 
more  or  less  velocity  modulation  may  Be  produced  depending  on  the 
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relaxation  constants  of  the  liquid.  Hov«ver  is  the  analysis 

took  velocity  modulation  asa  b-<sis,  this  did  not  ap  \y  to  the  theory, 

but  would  •’■'Vu  affcctc  \  the  experimental  results  somewhat . 

In  the  re  [d  on  of  predominant  bunching  and  disc  formation, 

the  restriction  1 5  low  velocity  laminar  jets  is  not  necessary- 

The  results  for  w'.tcr  at  95  and  145  fps  are  shown  in  Plates  IV. 

ani  V.  (taken  under  2(jO/;*  s  of  vibratj.cn  at  4,500)  cos). 

1,800 

Sore  a  line  ,-f  the  li  mi  !  is  caused  by  aerodynamic  fragmentation, 

but  the  definite,  mass  concentration  spaced  at  wave-length 

'istances  along  the  jot,  is  still  very  evident. 

It  is  of  interest  t''  note  that  single  discs  have 

(Q) 

been  observed  nreviously  on  liquid  jets.  They  were  induced  by  the 

application  of  a  steo  change  in  velocity  at  the  orifice.  The 
analysis  given  was  rather  crude  ani  wouj.i  not  cover  the  low 
modulation  velocity  region.  However  the  evidence  for  velocity 
modulation  at  th_  orifice  being  the  cause  of  the  phenomena 
renorted  here,  must  now  be  recorded  as  complete. 

Prediction  of  the  magaitude  (and  phase)  of  velocity 
modulation  produced  by  a  given  vioration  is  a  necessary  next  step. 

Prom  the  vibration  point  of  view  it  will  render  the  analysis  mor« 
complete.  Further  investigation  of  the  effect  of  vibration  on 
the  injector  cavity  hydrodynamics  is  planned. 
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PROPERTIES  OF  TAYLOR  VORTICES  OF  LARGS  AMPLITUDE 

This  paper  is  concerned  with  the  secondary  vortei 
flow  of  fluid  between  two  concentric  cylinders;  the  inner 
cylinder  rotates  at  a  constant  speed  while  the  outer  is 
kept  fixed.-  Taylor  (1923)  has  shown  that  as  the  3peed  of 
the  inner  cylinder  is  increased  beyond  a  certain  value  the 
motion  changes  from  a  purely  tangential  flow  to  one  in 
which  there  exists  a  superposed  secondary  vortex  flow.  The 
vertices  in  this  flow  are  in  the  shape  of  toroids  whose 
axis  is  the  axis  of  the  cylinders.  Taylor  predicted  the 
conditions  under  which  this  transition  takes  place;  he 
■also  gave  the  form  of  the  infinitesimal  vortices.  The 
extension  to  rates  of  rotation  above  the  critical  value,  when 
the  amplitude  of  the  secondary  flow  is  an  appreciable  fraction 
on  the  main  flow,  was  made  by  Stuart  (1958) . 

Taylor  confirmed  experimentally  his  theoretical 
value  for  the  critical  speed;  his  value  has  also  been 
verified  by  the  experiments  of  '•v'endt  (1933)  and  Donnelly  and 
Simon  (I960)  who  also  confirmed  Stuart’s  analysis.  Donnelly 
and  Simon  show,  however,  that  Stuart’ s  analysis  holds  only 
for  a  limited  range  of  speeds  above  the  critical.  It  is  the 
aim  of  this  paper  to  give  an  approximate  analysis  of  the 
flow  when  the  rates  of  rotation  are  so  large  that  Stuart's 
analysis  is  no  longer  valid. 

The  reason  for  the  occurence  of  these  vortices  lies 
in  the  centrifugal  force  field  set  up  by  the  primary  flow. 
When  the  notion  is  purely  tangential,  the  speed  of  the  fluid 
decreases  linearly  outwards  from  the  peripheral  speed  of 


/  A 


the  inner  cylinder  to  rest  on  the  outer  cylinder.  There 
exists,  therefore,  a  tendency  for  a  fluid  particle  to 
migrate  from  the  inner  cylinder  to  the  outer.  This  tendency 
will  be  realised  when  the  centrifugal  force  is  sufficient 
to  overcome  the  viscous  forces  in  the  fluid.  The  ratio  of 
these  forces  is  given  essentially  by  the  non-dimensional 


Taylor  number  T,  defined  by 

V  —  4?^ 


fT 


'  R  V  V  ' 

where  is  the  radius  of  the  inner  cylinder,  is  the 

gap  between  the  cylinders  and~U"  is  the  peripheral  speed 
of  the  inner  cylinder,  and  V  is  the  kinematic  viscosity. 

Taylor  showed  that  vortices  were  set  up  when  the 
Taylor  mmber  was  greater  than  about  1700.  When  T  is  greater 
than  1700  a  range  of  vortex  shapes  will  be  amplified;  the 
form  which  has  the  greatest  amplification  has  an  approximately 
square  cross-section,  this  form  ultimately  predominates. 

In  this  paper  it  is  shown  by  an  approximate  method 
that  the  drag  on  the  inner  cylinder  is  increased  by  these 

vortices  in  the  ratio  » / 

m  % 

7 

when  T>/  lo4  .  This  result  agrees  well  with  the 
experimental  results  of  Taylor  and  Wendt  quoted  above,  provided 
the  Reynolds  number  ■Ui/y  is  less  than  about  5,000, 
above  which  value  the  flow  becomes  turbulent. 


References 


Donnelly, R.J.  and  Simon  N.J.  (I960)  J.  Fluid  Mech.  7,401 
Stuart,  J ,?.  (1953)  J.  Fluid  Mech.  4,  1. 

Taylor,  G.I.  (1923)  Phil.  Trans.  A. 223,  289. 

Wendt,  F.  (1933)  Ingen.  -Arch.  4,577. 


f 


? 

i 


f 


* 


f 


-  41  - 


Part  3  (b) 


TITLE: 


"The  Tangential  Mode  and  Rocket  Engine  Burn-Out" 

Authors:  DR.  L.J.  CRIME 

®.  S.  BIRCH 
OR.  P.D.  McCORKACK. 


i 


i 

\ 


I! 

\) 


i 

» 


* 

l 

< 

? 

t 

« 

i 


f 

» 

t 

* 

f ; 

i ' 


Nomenclature ; 


^tvN  - 

^hVu»w,^  * 

03  - 

Ss 

*?  . 

Re  «  ^$/y 

t 

H'  * 

R  « 

F 

* 

X  - 


free  stream  velocity 

x  -  exponent  of  velocity  in  the  boundary  layer 

in  the  absence  of  vortices 

the  additional  velocity  components  due  to  the 

disturbance 

the  x  -  components  of  velocity  at  the  edges  of 
the  regions  of  abrupt  change 
vorticity  in  the  x  -  direction 
thickness  of  the  boundary  layer 
kinematic  viscosity 
«  Reynolds  Number 

thickness  of  the  region  of  abrupt  change  in  the  cell 

p3eudc  stream  function 

radius  of  curvature  of  the  wall 

skin  friction  ratio 

amplification  factor  for  snail  amplitude  values 
wave-length  in  the  z-direction 
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INTRODUCTION 

It  is  contended  here  that  the  impingement  of 
rotating  gas  on  the  wall  of  a  cylindrical  rocket  engine 
chamber,  as  will  occur  in  the  presence  of  tangential 
(although  this  involves  an  oscillating  flow)  and  spinning 
modes  of  instability,  will  result  in  the  three-dimensional 
type  of  boundary  layer  instability  first  designated  by 
Taylor  ^  ^  for  the  case  of  rotating  concentric  cylinders 

(3) 

and  later  by  G&rtler  for  gas  flow  over  a  concave  surface. 

The  formation  of  the  regular  and  well-ordered  pattern 
of  Taylor-Gortler  vortices,  with  axes  parallel  to  the  chamber 
wall,  will  be  snown  to  result  in  a  very  large,  albeit  oscil¬ 
lating,  increase  in  heat  flux  to  the  wall  over  that  due  to 
normal  convection. 

The  coordinate  system  used  by  Gortler  (see  "Boundary 
Layer  Theory"  by  H.  Schlichting  p.442)  will  be  adopted  in 
this  work, 

x  •*  coordinate  in  the  cylinder  wall  and  parallel  to  the 
main  gas  flow. 

y  a  coordinate  perpendicular  to  the  wall 
z  *  coordinate  perpendicular  to  the  x  direction. 

ll 

Gortler  analysed  the  time  build-up  of  the  vortices  by 
considering  a  small  three  dimensional  disturbance  superposed 
on  the  flow: 

u’  -  ~,(y)  cos  (JL'J  ) 
v’  »  v,(y)  cos  (JL,y)6^^ 
w‘  *  wf(y)  sin  (ot^ 

Instability,  or  vortex  formation,  will  be  initiated  when,  for 
a  given  wave  number  ,  the  characteristic  parameter  or 
Gortler  Number,  is  above  a  certain  critical  value: 

>n 


We  have  taken  the  experimental  results  quoted 
(4) 

by  H„C.  Kricg  Jnr.  for  a  15  inch  diameter  engine,  as 
a  basis  for  numerical  results. 

/  p\ 

Using  Schlichting' s  formula  ^  '  for  the  boundary 

layer  on  a  smooth  plate,  and  taking  a  free  stream  velocity 

-2 

of  2000  fps  and  a  maximum  available  time  of  10  sec,  a 

A 

mean  boundary  layer  thickness,  $  ,  of  about  3.5  x  10~J  ft. 


is  obtained. 

For  Krieg1 s  chamber, R  -  7.5  inches  and  so 
Gttrtler  dumber  ■  "HisJL  ,  »  46 

The  condition  for  instability  is  thus  satisfied. 


An  arbitrary  initial  disturbance  will  have  a 
range  of  Fourier  components  wnich  will  be  amplified.  The 
component  with  the  greatest  amplification  factor  corresponds 
to  a  cell  structure  whc3e  wave-length  is  about  *s<3).  Thi3 
cell  size  will  ultimately  predominate . 

In  the  case  considered  above,  j|J  for  this  type  of 
cell  is  about  40,000.  Thus,  the  initial  rate  of  build-up 
is  about  Kirchgassner  ^  has  proved  that  a 

finite  upper  bound  for  the  vortex  growth  exists,  this  being 
when  the  vortex  velocity  equals  the  free  stream  rotational 

velocity.  It  is  thus  certain  that  the  vortices  will  have 

-2 

reached  their  maximum  size  in  a  time  well  under  10  sec. 

The  ultimate  (steady-stat~'  rtex  s+^ncture  will  have  a 
wave-length  of  about  2  ^  . 

Analysis 

The  equations  of  motion  as  developed  by  GttrtXer, 
and  neglecting  terms  containing ,  are  as  follows: 


are  as  follows: 


c  v5i  +  ^  I  =  v  V  A* 


—  (1) 
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Uu,  _ii£<Vw  XVW1/Ul+jl„WVlV4^  (2) 
^  V£  +R^V  '  Y  . 

CJ  =  V1  —  (3) 

where  ^  is  a  pseudo  stream  function,  defined  by 


and  Cjis  the  component  of  vorticity  in  the  x  direction, 


to  = 


The  boundary  conditions  over  a  singLe  cell  (one  wave-length, 
two  vortices)  are  as  shown  in  Figure  1. 

Note:  all  the  quantities  are  periodic  in  the  z  direction, 

this  means  in  particular  that  A  ,  V( ,  Cj  have  the  same 
values  at  the  outside  edges  of  the  cell  respectively. 
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9y  an  order  of  magnitude  procedure  it  will  be 
shewn  that  the  rate  of  change  of  the  various  parameters, 

Cj  etc,,  cannot  be  of  the  order  ^/c  etc.,  over  ail- 

♦  ° 

the  cell. 

CJ  . 

<i  Suppose  the  variation  is  etc. 

Consider  equation  (2):  Term  by  term  this  leads  to  Table  1 


Table  1 


>—■  - . . . -  •  ■■■  - - 

•v  bV. 

•v  iii 

1  S 

J  JU . 

J  •  T 

y  S1 

IMif 

( 

V,  R 

57  T- 

Note:  Co  V/g 

Consider  equation  (1) :  Term  by  term  this  leady  to  Table  11. 


Table  11 


WA Mo  ^  ^Jix, 

— 

v-Vg> 

Hence  \T  *  fT~  ,  where 

Kxu 

Putting  this  value  for  *15",  ,  into  the  expressions  in  Table  1, 
the  bottom  row  new  becomes,  as  shown  in  Table  111. 


Table  111 


* 

I 

i 

y 


i 


s 

* 


I 

I 


: 
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The  first  and  third  terms  are  vanishingly  small 
(10‘5) ,  implying  that  the  boundary  conditions  cannot  be 
satisfied.  We  must  infer  then,  that  CJ  ,  jUL|  ,  etc.  only 
change  slovdy  over  most  of  the  cell,  there  being  a  thin 
region  where  very  abrupt  parameter  chancres  occur. 

Let  Its*  width  be  designated  by  the  symbol  £.  .  The  cell 
(see  Pig.  2)  now  consists  of  two  regions,  'slow*  and  ’abrupt* 
respectively. 

In  the  slew  region  the  viscous  terms  in  the 
equations  can  be  neglected.  New  "U0  » /&fl(y) ,  and  so 
equation  (1)  can  be  written  as, 


"t  n  =0orJlcnS^ 

where “jULo-J-jU,  ♦  This  implies  that  AX. 

so  that  in  the  ’slow*  region,  ( 0-+ jU  ,  )  is  a  constant  along 

a  stream-line. 

Derivation  of  £  in  terms  of  S  : 

In  the  ’abrupt*  region  we  have  from  equation  (1), 

«v .  2*  -  V-1*- 
1  S  £' 

Note:  the  variation  of  parameters  as  l/£  or  i/a  can  be 
assigned  here. 

Now  ,  » 

CO  -  ? 


and  so 


JL  ju.S* 


(6) 
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Since  it  has  been  oroved  that  ii  alone  a  stream- 

*  -  - — —  «_»  —  —  — 

line  in  the  inviscid  region  is  a  constant,  let  have  the 
constant  values  ^  streamlines  just  outside 

the  region  of  abrupt  change,  in  the  two  halves  of  the  cell 
res  actively. 


The  x-component  of  momentum  transferred  from  the 
main  stream  to  the  cell  must  bo  equal  to  the  skin  friction 
at  the  wall.  Hence, 

w 

•  •  i_  4-  -U  ^  ^  Uco 

As  A.tjto  near  the  wall,  jx  5  has  the  values  JU.  ^ 

respectively,  near  the  wall. 

Then, 


- (9) 


"\Xa3 


In  the  example  taken  here,  we  find  that  equation  9  leads  to 

a  value  of  F  of  about  10.  That  is,  a  10  fold  increase  in 

heat  flux  to  the  wall  is  predicted. 

Comparison  with  Published  Experimental  Results: 

The  only  extensive  set  of  data  for  heat  transfer 

thro1  a  fluid  between  concentric  rotating  cylinders  available, 

(7) 

is  that  of  Bjorklund  and  Kays  .  This  data  has  been  cor¬ 
related  by  Nissan  &  Haas  ^  .  —  - 
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Our  example,  with  a  critical  P  parameter  value 
of  12  (as  per  Nissan  and  Haas),  gives  a  P/P  of  about  200, 
and  thus 

Nusselt 

______  *  11 

(Nusselt) 

cond. 

(o) 

Similarly,  from  the  work  of  Donnelly,  Simpson  &  Batchelor  '  ' 
on  transfer  of  mcnentun  between  concentric  rotating  cylinders 
separated  by  a  low  viscosity  fluid,  the  fractional  increase 
in  mcmentum  transferred  (due  to  Taylor-Gftrtler  vortex  formation) 
is  predicted  to  be  of  the  order  of  20. 

We  understand  that  observations  on  heat  transfer 
to  the  rocket  engine  wall  in  the  pressence  of  tangential 
instability  indicate  a  fractional  increase  in  the  region 
of  ten. 

We  may  conclude  therefore,  that  experimental  evidence 
supports  the  contention  that  the  tangential  mode  of  instability 
in  a  rocket  engine  can  initiate  this  type  of  disturbance  in 
the  boundary  layer  on  the  wall,  and  lead  to  an  independent 
-and  significant  increase  in  heat  transfer. 
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Part  U. 

Future  Work 

An  experimental  study  of  a  gaseous  diffusion 
flame  will  be  carried  out  in  which  the  fuel  and  oxidant 
will  be  injected  coaxially  and  the  injector  plate  is 
vibrated-  transversely.  Schliercn  and  spectrometric 
techniques  will  be  used,  and  luminosity  and  pressure 
measurements  made  in  the  combustion  chamber,  to  3tudy  the 
flame  effects  induced  by  fuel  flow  modulation.  The  use 
of  gas  jets  rather  than  liquid,  eliminates  complications 
due  to  liquid  jet  break-up  and  evaporation. 

A  theoretical  analysis  cf  the  f lowing/ reacting 
system  in  the  burner  will  also  be  carried  out,  in  which 
the  gas  flew  is  modulating  with  finite  anplitude  and  at 
the  applied  vibration  frequency.  The  initial  analysis 
of  a  cold  jet  with  flow  modulation,  injected  into  still 
air,  has  almost  been  completed.  The  next  step  will  be 
to  consider  coaxial  jets  and  then  to  include  reaction  and 


the  chemical  kinetics. 


